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NUMERICAL SEMIGROUPS WITH A GIVEN SET OE PSEUDO-EROBENIUS NUMBERS 


M. DELGADO, P. A. GARCIA-SAnCHEZ, AND A. M. ROBLES-PEREZ 


Abstract. The pseudo-Frobenius numbers of a numerical semigroup are those gaps of the numerical semi¬ 
group that are maximal for the partial order induced by the semigroup. We present a procedure to detect if a 
given set of integers is the set of pseudo-Erobenius numbers of a numerical semigroup and, if so, to compute 
the set of all numerical semigroups having this set as set of pseudo-Frobenius numbers. 


1. Introduction 

Let S be a numerical semigroup, that is, a cofinite submonoid of (IN,-!-), where IN stands for the set of 
nonnegative integers. 

An integer x is said to be the Frobenius number of S (respectively, a pseudo-Frobenius number of S) if 
X g S and x -i- s € S, for all s € IN \ {0} (respectively, for all s € S \ {0}). 

Given a positive integer /, there exist numerical semigroups whose Frobenius number is /. One 
example of such a semigroup is the semigroup {0,/ + 1, —>} containing 0 and all the integers greater than 
/. There are several algorithms to compute all the numerical semigroups with a given Frobenius number 
(the fastest we know is based on [9]). 

We denote by F(S) the Frobenius number of S and by PF(S) the set of pseudo-Frobenius numbers of 
S. The cardinality of PF(S) is said to be the type of S and is denoted by t(S). 

A positive integer that does not belong to S is said to be a gap of S and an element of S that is not 
greater than F(S)-i-1 is said to be a small element of S. To denote the set IN\ S of gaps of S we use gaps(S) 
and to denote the set of small elements of S we use smalls(S). Since a set of gaps must contain the 
divisors of all its members and a set of small elements must contain all multiples of its members (up 
to its maximum), it is clear that there are sets of positive integers that can not be the set of gaps or the 
set of small elements of a numerical semigroup. The set of gaps, as well as the set of small elements, 
completely determines the semigroup. Observe that when some elements or some gaps are known, 
others may be forced. For instance, a gap forces all its divisors to be gaps. 

Let « be a positive integer and let PF = {gi,g 2 >--->gn-i>Sn} be a set of positive integers. Denote by 5(PF) 
the set of numerical semigroups whose set of pseudo-Frobenius numbers is PF. When « > 1, 5(PF) may 
clearly be empty. Moreover, when non-empty, it is finite. In fact, 5(PF) consists of semigroups whose 
Frobenius numbers is the maximum of PF. Some questions arise naturally. Among them, we can consider 
the following. 

Question 1. Find conditions on the set PF that ensure that 5(PF) ^ 0. 

Question 2. Find an algorithm to compute 5(PF). 

Both questions have been solved for the case that the set PF consists of a single element (which must be 
the Frobenius number of a numerical semigroup; symmetric numerical semigroups) or when PF consists 
on an even positive integer / and //2 (pseudo-symmetric numerical semigroups), see [1]. Moreover, 
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Question 1 was solved by Robles-Perez and Rosales [7] in the case where PF consists of 2 elements (not 
necessarily of the form {/,//2}). 

The set 5(PF) can be computed by filtering those semigroups that have PF as set of pseudo-Frobenius 
numbers, from the numerical semigroups whose Frobenius number is the maximum of PF (cf. Exam¬ 
ple 12). Due, in part, to the possibly huge number of semigroups with a given Frobenius number, this is 
a rather slow procedure and we consider it far from being a satisfactory answer to Question 2. 

Irreducible numerical semigroups with odd Frobenius number correspond with symmetric numer¬ 
ical semigroups, and those with even Frobenius number with pseudo-symmetric (see for instance [8, 
Chapter 3]). Bresinsky proved in [2] that symmetric numerical semigroups with embedding dimension 
four have minimal presentations of cardinality 5 or 3 (complete intersections). Symmetry of a numeri¬ 
cal semigroup S translates to having {F(S)} as set of pseudo-Frobenius numbers. Eater Komeda, [6], was 
able to prove the same result for pseudo-symmetric numerical semigroups (though he used different ter¬ 
minology for this property; in this setting 3 does not occur since pseudo-symmetric are never complete 
intersections). A numerical semigroup S is pseudo-symmetric if its set of pseudo-Frobenius numbers 
is {F(S),F(S)/2}. It should be interesting to see the relationship with the type and the cardinality of a 
minimal presentation, and thus having tools to find semigroups with given sets of pseudo-Frobenius 
numbers becomes helpful. Watanabe and his students Nari and Numata are making some progress in 
the study of this relationship. 

1.1. Contents. We present two different procedures to determine the set of all numerical semigroups 
with a given set of pseudo-Frobenius numbers. One exploits the idea of irreducible numerical semi¬ 
group. From each irreducible numerical semigroup we start removing minimal generators with certain 
properties to build a tree whose leafs are the semigroups we are looking for. The other approach is 
based on determining the elements and gaps of any numerical semigroup with the given set of pseudo- 
Frobenius numbers, obtaining in this way a list of "free" integers. We then construct a binary tree in 
which branches correspond to assuming that these integers are either gaps or elements. 

We start this work with some generalities and basic or well known results and connections with the 
GAP [5] package numericalsgps [4] (Sections 2 and 3). Then we describe a procedure to compute forced 
integers (Sections 4 and 5). As computing forced integers is fast and leads frequently to the conclusion 
that there exists no semigroup fulfilling the condition of having the given set as set of pseudo-Frobenius 
numbers, this approach benefits from the work done in the preceding sections. The following section of 
the paper is devoted to the above-mentioned approach of constructing a tree with nodes lists of integers, 
which turns out to be faster most of the times than the one based on irreducible numerical semigroups. 
Nevertheless, besides being useful to compare results, the method using irreducible numerical semi¬ 
groups is of theoretical importance so we decided to keep it, and it is given in Appendix A. The next sec¬ 
tion describes an algorithm that (increasing the number of attempts, if necessary) returns one numerical 
semigroup with the given set of pseudo-Frobenius numbers, when such a semigroup exists. Although 
it may be necessary to increase the number of attempts, usually does not require the computation of all 
the semigroups fulfilling the condition and is reasonably fast in practice (Section 8). 

We give pseudo-code for the various algorithms that are used in the approach of lists of free integers. 
The pseudo-code given has an immediate translation to the GAP language, which is the programming 
language we used to implement these algorithms. We observe that it is a high level programming lan¬ 
guage. Note that we take advantage of the existence of the GAP package numericalsgps for computing 
with numerical semigroups. The names used for the functions described here are slightly different from 
the ones used in the package, since there longer names are required to be accurate with the variable 
names policy in GAP. Essentially, it will be a matter of adding a suffix to the names used for the package. 

Many examples are given throughout the paper, some to illustrate the methods proposed, while others 
are included to motivate the options followed. In some of the examples we show the output obtained in a 
GAP session. These usually correspond to examples that are not suitable to a full computation just using 
a pencil and some paper; furthermore, indicative running times, as given by GAP, are shown (mainly in 
Section 7). 

A new version of the numericalsgps package, including implementations of the algorithms developed 
in the present work, is released at the same time this work is made public. The implementations can 
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be checked (the software in question is open source; see the links in the references) and can be used for 
testing examples. 


2. Generalities and basic results 

Throughout the paper we will consider often a set PF = {gi,g 2 ^---^gn-i^gn} of positive integers. We will 
usually require it to be ordered, that is, we will assume that gi <g 2 < ■ ■ ■ < g„-i < gn- For convenience, we 
write PF = < g 2 < • • • < g„-i < gn} in this case. 

We denote by frob the maximum of PF and by type the cardinality of PF. Note that if S € 5(PF), then 
frob = g„ = F(S) and type = n = t(S). 

2.1. Forced integers. We say that an integer ^ is a gap forced by PF or a P¥-forced gap if ^ is a gap of 
all the numerical semigroups having PF as set of pseudo-Frobenius numbers. In particular, if there 
is no semigroup with PF as set of pseudo-Frobenius numbers, then every nonnegative integer is a gap 
forced by PF. We use the notation QT{PV) to denote the set of PF-forced gaps. In symbols: QT{PV) = 

ns65(PF)gaps(S). 

In a similar way, we say that an integer £ is an element forced by PF or a PP-forced element if £ is an 
element of all semigroups in 5(PF). We use the notation £PF{PP) to denote the set of (small) PF-forced 
elements. In symbols: £T(PP) = ns 65 (PF)S^alls(S). Note also that if 5(PF) = 0, then £T(PP) = IN. 

The union of the PF-forced gaps and PF-forced elements is designated by PP-forced integers. 

The following is a simple, but crucial observation. 

Proposition 3. 5(PF) ^ 0 if and only ifQF{PP) n £F{PP) = 0. 

Proof If 5(PF) = 0, then all nonnegative integers are at the same time gaps and elements forced by PF. 
Conversely, assume that 5(PF) ^ 0 and let S e 5(PF). Then QT(PP)r\£T(PP) c gaps(S)nsmalls(S) = 0. □ 

Frequently the prefix PF is understood and we will abbreviate by saying ]\isi forced gap, forced element 
or forced integer. 

Let G and E be, respectively, sets of forced gaps and forced elements. The elements v e {!,...,frob} 
that do not belong to G U £ are said to he free integers for (G,E). When the pair (G,E) is understood, we 
simply call free integer to a free integer for (G, E). 

2.2. Well known results. The partial order <5 induced by the numerical semigroup S on the integers 
is defined as follows: x <5 y if and only if y - x € S. The following result is well known and will be used 
several times throughout this paper. 

Lemma 4. [8, Lemma 2.19] Let S be a numerical semigroup. Then 

(i) PF(S) = Maximals<j(Z\S), 

(ii) X € Z \ S if and only if f -x e S for some f € PF(S). 

It is well known that the type of a numerical semigroup S is upper bounded by the least positive 
element belonging to S, which is known as the multiplicity of S and it is denoted by m(S). 

Lemma 5. [8, Corollary 2.23] Let S be a numerical semigroup. Then m(S) > t(S) -F 1. 

2.3. Initializing the set of forced gaps. The maximality of the pseudo-Frobenius numbers of S with 
respect to <5 means that they are incomparable with respect to this ordering. In particular, the difference 
of any two distinct pseudo-Frobenius numbers does not belong to S, that is, 

(1) [gi-gj I hi €{l,...,t(S)],i >y]Cgaps(S). 

This is the underlying idea of the next result. 

Lemma 6. Let S be a numerical semigroup and suppose that PF(S) = {gi < gi ^ ^ gn-i < gnl^ with n> 1. 

Let i € {2, ...,t(S)] and g e (PF(S)) with g < gi. Then gi-g ^ gaps(S). 
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Proof. Assume that g = H-for some fc € IN. We proceed by induction on k. The case fc = 1 is given 

by (1). Assume that the result holds for k-l and let us prove it for k. If gi-g ^ S, then + {gj -g) e S 

by definition of pseudo-Frobenius number. It follows that gj - (gj^ + -J € S, contradicting the 

induction hypothesis. □ 

Remark 7. Lemma 5 implies that {x € IN | 1 < x < t(S)} c gaps(S). Hence c ^jr(PF). 

As the pseudo-Frobenius numbers of S are gaps of S and any positive divisor of a gap must be a gap 
also, we conclude that the set of divisors of 

PF(S) U {x € N 11 < X < t(S)} U [g,-g\ie {2,...,t(S)},g e {PP(S)fg, > g} 

consists entirely of gaps of S. 

Consider the set 

PFU {X€ N I 1 <x < «} U {g; -g I z € {2,...,n],ge{PV),gi >g} 

and denote by sfg(PF) the set of its divisors (as we are only considering positive divisors, in what follows 
we will not include this adjective). If S is a numerical semigroup such that PF(S) = PF, we deduce that 
sfg(PF) c gaps(S). We have proved the following result for the case where there is a numerical semigroup 
S such that PF(S) = PF. If no such semigroup exists, then GJ^{PF) = IN and the result trivially holds. 

Corollary 8. Let PF be a set of positive integers. Then sfg(PF) c Q7F{PV). 

We use the terminology starting forced gap for PF to designate any element of sfg(PF), since sfg(PF) is 
the set we start with when we are looking for forced gaps. In Subsection 4.1 we provide pseudo-code for 
a function to compute starting forced gaps. 

2.4. Initial necessary conditions. Let « > 1 be an integer and let PF = {gi < • • • < gn} be a set of positive 
integers. 

Lemma 9. Let S be a numerical semigroup such that PF(S) = PF. Let i e {2,...,«} and g € (PF) \ {0} with 
g < gj. Then there exists A: € {!,...,«} such that gj^ - {gj -g) e S. 

Proof. Lemma 6 assures that gj- g€S. The conclusion follows from Lemma 4. □ 

By choosing z = « in the above result, there exists k^ n such that gfc - (gn -gi) > 0 and gi^ - {gn - gi) ^ 
{!,...,type}. (Note that k = n would imply gi € S, which is impossible.) But then g„_i -(gn“gi) ^ 0^ since 
gn-i ^ gk for all A: € {1,..., zz - Ij. We have thus proved the following corollary. 

Corollary 10. Let S be a numerical semigroup such that PF(S) = PF. Then gi > g„ -g„-i. 

The computational cost of testing the condition gi > g„ - g„-i obtained in Corollary 10 is negligible 
and should be made before calling any procedure to compute 5(PF), avoiding in many cases an extra 
effort that would lead to the empty set. 

Other conditions of low computational cost would also be useful. Since g„_i-{g„-gi) > 0, one could be 
tempted to ask whether replacing gi by gj one must have g„_i - {gn-gi) ^ {L---Aypej (since {1,.. .,type} 
consists of gaps). The following example serves to rule out this one that seems to be a natural attempt. 

Example 11. Let S = (8,9,10,11,13,14,15). One can check easily that PF(S) = {5,6,7,12}. For PF = 
{5,6,7,12}, we have g„_i -(g„-g 2 ) = 7-(12-6) = 1 € {!,...,4} = {1,...,type}. 

Later, in Subsection 4.3, we give an extra condition, which is based on forced integers. 

3. Connections with GAP and the package numericalsgps 

When developing the present algorithms we had one implementation in mind. As already referred, 
it was to be made in the GAP programming language, benefiting also of the package numericalsgps. In 
fact the implementation accompanied the development of the algorithms, with mutual benefits. 

To understand the paper no previous familiarity with GAP is assumed. We give some examples to get 
the familiarity needed to fully understand the pseudo-code. Many of the terms used for the pseudo¬ 
code presented are borrowed from the GAP programming language, which is a high level one. The terms 
Union, Length, Difference are used with the meaning they have in GAP, which is clear. 
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Like what is done in GAP, for an unary function that applied to an argument arg returns an expression 
expr, we write 


arg expr. 

This shorthand for writing a function is used in our pseudo-code. Next we briefly explain the meaning 
of some of the terms used. For a not so brief explanation, complete definitions and plenty of examples, 
GAP's manual should be consulted. 

Here list is a list (of integers or of numerical semigroups) and func is an unary function (that applies 
to the objects in list). 

• List{list, func) returns a new list new of the same length as list such that new[i] = func(list[i]). 

• Filtered(/zst,/u«c) returns a new list that contains those elements of list for which func re¬ 
turns true. The order of the elements in the result is the same as the order of the corresponding 
elements of this list. 

• One can use Set{list) (which is a synonym of SSortedList ("strictly sorted list")) to get a list that 
is duplicate free and sorted (for some order). 

• M6Set{list,el) adds the element el to the set list. 

• F icst{list,func) returns the first element of list for which the unary function func returns true. 
If func returns false for all elements of list, then First returns fail. 

• ForAll(/zst,/u«c) tests whether func returns true for all elements in list. 

• ForAny(/zst,/u«c) tests whether func returns true for at least one element in list. 

• IsRange(/zst) detects if the argument is an interval of integers. 

We also use the following abbreviations. Let X be a list of integers and P a list of positive integers. 

• Posint (X), an abbreviation of Filtered(X, IsPosInt), returns the positive integers of X. 

• Divisors (P), an abbreviation of Union (List (P ,DivisorsInt)), returns the divisors of the ele¬ 
ments of P: it first computes the lists of divisors for each element in P, and then takes the union 
of all of them. 

The numericalsgps package also influences our pseudo-code. In some cases we use directly the names 
of the available functions, but in some cases we use shorthands, which are intended to turn the pseudo¬ 
code more readable. 

One of the functions we use is RepresentsGapsOf NumericalSemigroup, which, for a given set X of 
positive integers returns true in case there exists a numerical semigroup S such that gaps(S) = X, and 
returns/fl/se otherwise. 

The functions to produce numerical semigroups used here are NumericalSemigroup, which is used 
when generators are given, and NumericalSemigroupByGaps, which, for an input X, returns a numerical 
semigroup whose set of gaps is X, when such a semigroup exists. 

As a shorthand for SmallElements(S), which gives the elements of the numerical semigroup S that 
are nor greater than F(S) -F 1, we use simply smalls, which agrees with the notation already introduced. 
Closure(elts,frob) is a shorthand for 

NumericalSemigroup(Union(elts,[frob+1..frob+Minimum(elts)])) 

which gives the least numerical semigroup containing the set elts of positive integers and having the 
largest possible Frobenius number not greater than frob. Note that the minimum of elts is greater than or 
equal to the multiplicity of the semigroup. Therefore, the union considered ensures that the semigroup 
contains all integers that are greater than the number frob given (aiming to be an upper bound for the 
Frobenius number of the semigroup constructed, although in cases where the elements elts define a 
numerical semigroup, it may have smaller Frobenius number). 

As in GAP, comments start with the character #. 

As an example on how we will present function/algorithms in this manuscript, we write the function 
Closure. 
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Closure(elts,frob) # the least numerical semigroup containing elts and 

# whose Frobenlus number Is not greater than frob 
return NumerlcalSemlgroup (Union (elts, [frob +l..frob +Mlnlmum (Posltlvelnt (elts))])); 

Function Closure 

The following example is just a GAP session that is intended to illustrate how to compute the set 
5(PF) by filtering those semigroups that have PF as set of pseudo-Frobenius numbers, from the numer¬ 
ical semigroups whose Frobenius number is the maximum of PF. This process was mentioned in the 
introduction. 

Throughout the examples, we use a simple way provided by GAP to give a rough idea of the time 
spent: time is a global variable that stores the time in milliseconds the last command took. 

Example 12. We illustrate how to compute the set of numerical semigroups having {19,29} as set of 
pseudo-Frobenius numbers. Double semicolon in GAP inhibits the output. 

gap> pf := [19,29];; 

gap> nsf29 := NumerlcalSemlgroupsWlthFrobenlusNumber(29);;tlme; 

31372 

gap> Length(nsf29); 

34903 

gap> nspf1929 := Filtered(nsf29, s -> PseudoFrobenlusOfNumerlcalSemlgroup(s) = pf);; 
gap> time; 

2540 

gap> Set(nspf1929,MlnlmalGeneratlngSystem); 

[ [ 3, 22, 32 ], [ 6, 9, 16, 26 ], [ 7, 9, 17 ], [ 8, 9, 14 ], [ 8, 9, 15, 22, 28 ], 

[ 9, 12, 13, 14 ], [ 9, 12, 13, 15, 23 ], [ 9, 12, 14, 16, 22 ], 

[ 9, 12, 15, 16, 22, 23, 26 ], [ 9, 13, 14, 17, 21, 24, 25 ], 

[ 9, 13, 15, 17, 21, 23, 25 ], [ 9, 14, 16, 17, 21, 22, 24 ], 

[ 9, 15, 16, 17, 21, 22, 23, 28 ] ] 

Prior to the obtention of the procedures that are the object of study of the present paper, we got the 
necessary insight through detailed analysis of many examples with PF consisting of small numbers (less 
than 29, say). Example 12 illustrates that this can be easily done by using the numericalsgps package. 
When PF consists of small numbers the time spent is acceptable. 

The implementation of our algorithm (available in Version 0.99 of the numericalsgps package) is 
much faster. It is even faster than just filtering among the numerical semigroups with a given set of 
pseudo-Frobenius numbers. 

Example 13. To illustrate, we continue the GAP session started in Example 12. 

gap> new := NumerlcalSemlgroupsWlthPseudoFrobenlusNumbers(pf);;tlme; 

29 

gap> Set(new,MlnlmalGeneratlngSystem)=Set(nspf1929,MlnlmalGeneratlngSystem); 
true 

Visualization of images obtained by using the intpic [3] GAP package have also helped us to improve 
our understanding of the problem and to get the necessary intuition. 

The production of Figure 1 takes less than two minutes in our laptop (using numericalsgps and int¬ 
pic packages). Taking into account that there are 34903 numerical semigroups with Frobenius num¬ 
ber 29, one can consider the performance of the implementations satisfactory. In order to explain the 
meaning of the colors used in this picture to highlight some integers, we anticipate the results shown 
in Example 23. For PF = {19,29], we have that {1,2,4,5,10,11,19,20] consists of forced gaps and that 
{0,9,18,24,25,27,28,30] consists of forced elements. To the elements in each of these sets, as well as the 
ones in the sets of minimal generators, is assigned one color {red corresponds to pseudo-Frobenius num¬ 
bers, blue to minimal generators, green to elements, cyan to forced gaps, and magenta to forced elements; 
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Figure 1. The numerical semigroups with pseudo-Frobenius numbers {19,29}. 


in a black and white visualization of this manuscript this will correspond with different gray tonalities). 
For integers that belong to more that one set, gradient colors are assigned. 

4. Integers forced by pseudo-Frobenius numbers - some preliminary procedures 

In this section we give pseudo-code for several functions implementing results given in Section 2 and 
some others needed later in Section 5. 

Once more, PF is a fixed set {gi < g2< ■ < gn-i < of positive integers; frob stands for g„ and type 

stands for n. 

4.1. Forced gaps. The function StartingForcedGaps returns the integers considered in Subsection 2.3, 
which we called starting forced gaps. Recall (Corollary 8 ) that these have to be gaps of all the numerical 
semigroups having PF as set of pseudo-Frobenius numbers. The justifications for Lines 1 and 2 are given 
by Remark 7 and Lemma 6 , respectively. Line 3 is justified by Lemma 4(ii): when it is detected a gap 
that had to be an element (forced by exclusion), there exists a contradiction and the function returns fail. 
This proves the following proposition. 

Proposition 14. Let OUT be the output of the function StartingForcedGaps for the input PF. If OUT is not 
fail, then it consists of gaps of any numerical semigroup S € 5(PF). 


StartingF orcedGaps(PF) 

1 /^ := Union([l..type],PF); # uses the fact m(S) > t(S) -I- 1 (Remark 7) 

closures := List([l..type - l],i ^ smalls(Closure({PF[l],PF[2],...,PF[i]|,PF[i -h 1]))); 

diffs := 0; 

2 for i € [2..type] do 

^ := Union(dij5^s,PositiveInt(PF[i] -c/osMres[i - 1])); # see Lemma 6 

fg := Divisors(Union(^,dij5^s)); 

# Detect possible contradiction given by Lemma 4(ii) 

3 for X € Difference(/^, PF)) do 

if PositiveInt(Difference(PF-x,/^)) = 0 then 
L return fail 

return fg; 

Function StartingForcedGaps 


Example 15. Let PF = {16,29}. One can easily check that the set of starting forced gaps is {1,2,4,8,13,16,29}. 
The corresponding GAP session: 
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gap> pf := [16,29];; 

gap> g := StartingForcedGapsForPseudoFrobenius(pf); 

[ 1, 2, 4, 8, 13, 16, 29 ] 

The function FurtherForcedGaps is used to determine forced gaps when some gaps and some elements 
of a numerical semigroup are known. The justification for the fact that the output of FurtherForcedGaps 
consists of gaps (unless there is an element that also had to be a gap, in which case it returns/flz7) is the 
following: if f - e and e are elements of a semigroup then e + f - e = f belongs to the semigroup. In 
particular, if e is an element and / is a gap, then f - e is either negative or a gap. This proves the 
following proposition. 

Proposition 16. Let OUT be the output of the function FurtherForcedGaps for the input (fg, fe), with fg and 
fe consisting ofPV-forced gaps and P¥-forced elements, respectively. If OUT is fail, then 5(PF) = 0. Otherwise, 
OUT consists of gaps of any numerical semigroup S € 5(PF). 


FurtherForcedGaps (fg,fe) 

# fg and fe consist of gaps and elements. 


# respectively 

ng := Divisors(PositiveInt(Union(List(fg,g - 
1 if Intersection(fe, «g) = 0 then 


return Union(fg,«g); 


else 


^ return fail; 



Function FurtherForcedGaps 


4.2. Forced elements. We use two ways to get new forced elements. One of these ways makes use of 
Lemma 4(ii). We refer to the elements obtained in this way as elements forced by exclusion. Another 
way makes use of the following lemma, which tells us that small gaps force elements that are close to 
the maximum of PF. Sometimes we refer to them by using the more suggestive terminology big forced 
elements. 

Lemma 17. Let m be the multiplicity of a numerical semigroup S and let i be an integer such that I < i < m. 
Then either F(S) - i e S or F(S) - i € PF(S). 

Proof. It suffices to observe that, as i < m, one has that F{S)-i+m > F(S), and consequently F(S)-i + m e S. 
The result follows immediately from the definition of pseudo-Frobenius numbers. □ 

We observe that in the function FurtherForcedElements it is used that PF is precisely the set of 
pseudo-Frobenius numbers; otherwise there is no guarantee that the output consists of forced elements. 

Let us now prove the correctness of this function. Justification for the the result produced by the cycle 
starting in Line 2 is given by the direct implication of Lemma 4(ii). For the cycle starting in Line 3 is 
given by the reverse implication of the same lemma. (The integers known to be forced gaps are assumed 
to be gaps.) 

Justification for Line 1 follows from Lemma 17, since the m appearing there is smaller than or equal 
to the multiplicity. 

We have then the following proposition. 

Proposition 18. Let OUT be the output of the function FurtherForcedElements for the input (fg,fe), with fg 
and fe consisting of PP-forced gaps and PP-forced elements, respectively. If OUT is "fail”, then 5(PF) = 0. 
Otherwise, OUT consists of elements of any numerical semigroup S € 5(PF). 
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respectively 
a forced gap 

1 be := Dif ference(frob- - 1},PF); 

# Elements forced by exclusion 

ee := 0; 

2 for X € fg do 

/z7f := Filtered(PF,/^ not((/ -X € fg) or (f - x < 0))); 
if Fength(/z7t) = 1 and /z7f[l] - x ^ ee then 
^ AddSet(ee,/z7i[l]-x); 

candidates := Difference([l..frob- 1],Union(fg,Union(fe,re))); 

3 ior X e candidates do 

if PositiveInt(Dif ference(PF-x,fg)) = 0 then 
^ AddSet(ee,x); 

ne := \]nion{fe,ee,be); 

4 if Intersection(fg, «e) = 0 then 

return ne; 
else 

^ return fail; 

Function FurtherForcedElements 


FurtherForcedElements(fg,fe) # fg and fe consist of gaps and elements, 
# Big forced elements 

m := First(lntegers,« ^ > 0 and « £ fg); # least integer that is not 


4.3. A condition based on forced integers. When searching for forced integers, one should pay atten¬ 
tion to the existence of possible contradictions. 

Example 19. Fet PF = {4,9}. Taking divisors and the difference 9-4, one immediately sees that the set 
of starting forced gaps contains jl, 2,3,4,5,9}. But then 5 appears as a forced gap and as a (big) forced 
element. This is a contradiction which shows that {4,9} cannot be the set of pseudo-Frobenius numbers 
of a numerical semigroup. 

According to Proposition 3, the only that we need to take into account is that the set of forced gaps is 
disjoint from the set of forced elements. 

At the end of the functions (Fine 1 in function FurtherForcedGaps, and Fine 4 in function Further¬ 
ForcedElements), possible contradictions are detected. 

5. A PROCEDURE TO COMPUTE INTEGERS FORCED BY PSEUDO-FrOBENIUS NUMBERS 

Fet PF = {gi,g 2 >--->gn-i>gn} be a set of positive integers. The aim of this section is to give a procedure 
to compute some elements of QT{PV) and £T{PV). 

We start with a subsection that contains a procedure that makes use of the functions given in Section 4, 
which implement some well-known facts recalled in Section 2. 

We then devote a subsection to what we call admissible integers. This will yield another procedure 
(in Subsection 5.3), which, at the cost of increasing the execution time, may find more forced integers. 

The main algorithm of the present paper (Algorithm 3) would work as well by using only the quick 
version. Our experiments led us to consider the option of using the slower version once, and then use 
the quick version inside a recursive function that is called by the main algorithm. 

5.1. A quick procedure to compute forced integers. The correctness of Algorithm 1 (SimpleForced- 
Integers) follows from Propositions 16 and 18 which state the correction of the functions considered 
in Section 4. 

Theorem 20. Algorithm 1 correctly produces the claimed output. 
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Input : g, e, where g, e are sets of PF-forced gaps and PF-forced elements, respectively. 

Output: [fg, fe ], where fg 3 g is a set of PF-forced gaps and fe 3 e is a set of PF-forced elements, or 
fail, when some inconsistency is discovered 

fg := ShallowCopy(g); # used to store new forced elements and gaps 

fe := ShallowCopy(e); # without creating conflicts in the memory 

repeat 

changes := false; 

1 gaps := FurtherForcedGaps(fg,fe); 
if gaps = fail then 

return fail 
else 

if gaps ^ fg then 
changes := true; 

L fg-= gaps; 

2 elts := FurtherForcedElements(fg,fe); 
if elts = fail then 

return fail 
else 

if elts ^ fe then 
changes := true; 
fe := elts; 

until not changes; 
return [fg,fe]; 

Algorithm 1: SimpleForcedIntegers 


Note that the forced elements returned by Algorithm 1 (the list in the second component) are obtained 
by applying the function Closure to some set. Therefore we observe the following. 

Remark 21. The second component of the list returned by Algorithm 1 is the set of small elements of a 
numerical semigroup. 

Example 22. Let us see how the SimpleForcedIntegers works for PF = {16,29}. We already know (Ex¬ 
ample 15) that the set of starting forced gaps is sfg = jl, 2,4,8,13,16,29}. 

Let us now make a call to SimpleForcedIntegers with input sfg, [ ]. 

The first passage in Line 1 does not produce any new integer. 

The first passage in Line 2 produces 

jO, 3,6,9,12,15,18,21,24,25,27,28,30} 

as current forced elements. Observe that 3 is forced by exclusion (note that 3 = 16-13; also, 29 -13 = 16 
and 16 is a forced gap); 25 is also forced by exclusion (note that 16 - 25 < 0 and 29 - 25 = 4 is a forced 
gap). Also, 21 is forced by exclusion, but for now we do not need to worry with the multiples of 3, 
because these will appear when taking the closure. 

The second passage in Line 1 produces 

11,2,4,5,7,8,10,11,13,14,16,17,20,23,26,29} 

as current forced gaps. To check it, observe that 29 - 3 = 26, 29 - 6 = 23, 29 - 9 = 20, etc. are forced gaps. 
But then, 10 and 5 are (forced to be) gaps. 

No further forced elements appear. In fact, the union of the sets of forced gaps and of forced elements 
is jo,..., 30}. 

Therefore, all positive integers less than 29 are forced. One can check that the closure of the set of 
forced elements does not produce new forced elements, thus it is the set of small elements of a numerical 
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semigroup. Also, one can check that no forced gap outside {16,29} is a pseudo-Frobenius number, thus 
one may conclude that there exists exactly one numerical semigroup S such that PF(S) = {16,29}. 

This example illustrates that more than one passage through the repeat-until loop of the algorithm 
SimpleForcedlntegers may be needed. 

Example 23. In this example we present a GAP session, now for PF{19,29}. 
gap> pf := [19,29];; 

gap> sfg := StartingForcedGapsForPseudoFrobenius(pf); 

[ 1, 2, 5, 10, 19, 29 ] 

gap> SimpleForcedIntegersForPseudoFrobenius(sfg,[],pf); 

[ [ 1, 2, 4, 5, 10, 11, 19, 20, 29 ], [ 0, 9, 18, 24, 25, 27, 28, 30 ] ] 

Recall that the names used in our package are longer than the ones used in this manuscript, so for 
instance, SimpleForcedlntegersForPseudoFrobenius is the name we have used in our package for 
Algorithm 1. This example is related to Figure 1, where forced integers are highlighted. 

Example 24. Let us now apply the algorithm to PF = {15,20,27,35}. Again, we will use GAP to help us in 
doing the calculations (which can be easily confirmed by hand). 
gap> pf := [ 15, 20, 27, 35 ];; 

gap> sfg := StartingForcedGapsForPseudoFrobenius(pf); 

[ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 15, 20, 27, 35 ] 

We immediately get that {25,26,28,29,30,31,32,33,34} consists of forced big elements. And one can 
observe that 19 and 23 are forced by exclusion. This leads to the obtention of 35 - 19 = 16 as forced gap. 
No other forced elements are obtained, which agrees with the following continuation of the GAP session: 
gap> SimpleForcedIntegersForPseudoFrobenius(sfg,[],pf); 

[ [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 15, 16, 20, 27, 35 ], 

[ 0, 19, 23, 25, 26, 28, 29, 30, 31, 32, 33, 34, 36 ] ] 

5.2. Admissible integers. Let G and E respectively be sets of PF-forced gaps and PF-forced elements, 
and let v be a free integer for (G, £). We say that v is admissible for (G,E) if Algorithm 1 when applied 
to (G, £ U {v}) does not return fail. Otherwise, we say that v is non-admissible for {G,E). Thus, v is non- 
admissible implies that v can not be an element of any semigroup in 5(PF) and therefore is a gap of all 
semigroups in 5(PF), that is, is a forced gap. 

Lemma 25. Let G and E be sets of forced gaps and forced elements, respectively. Let v be free for (G,E). Ifv is 
non-admissible for {G,E), then v is a PF-forced gap. 

Observe that a semigroup generated by admissible elements for some pair (G,E) consists of admissible 
elements for (G,£). 

The function NonAdmissible, with input a pair of sets of forced gaps and forced elements, returns 
non-admissible integers, which, by Lemma 25, are new forced gaps. This function is called by Algo¬ 
rithm 2, a not so quick procedure to compute forced integers. 

Proposition 26. Let OUT be the output of the function NonAdmissible for the input (fg,fe), ivith fg and 
fe consisting of PF-forced gaps and PF-forced elements, respectively. Then OUT is a set of integers that are 
non-admissible for (fg,fe). 

Example 27. Let Pf = {11,22,23,25}. 
gap> pf := [ 11, 22, 23, 25 ];; 

gap> sfg := StartingForcedGapsForPseudoFrobenius(pf); 

[ 1, 2, 3, 4, 5, 6, 7, 11, 12, 14, 22, 23, 25 ] 

By using the function SimpleForcedlntegers one obtains the following. 
gap> SimpleForcedIntegersForPseudoFrobenius(sfg,[],pf); 

[ [ 1, 2, 3, 4, 5, 6, 7, 11, 12, 14, 22, 23, 25 ], 

[ 0, 18, 19, 20, 21, 24, 26 ] ] 
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NonAdmissible(fg,fe) # fg and fe consist of forced gaps and forced eiements, 

respectiveiy 

admissible := 0; 

totesf := Difference([l..frob],Union(fg,fe)); # start with the free integers 

while totest 0 do 
V := totest[l]; 

pnfce := SimpleForcedIntegers(fg,Union(fe, [v]), PF); 
if pnfce ^ fail then 

admissible := \]nion{admissible,pnfce[2]); 
totest := Di££ecence{totest,admissible); 

else 

^ totest := Dif f erence{totest, [v]); 

return Difference([l..frob], fldmzsszfo/e); 

Function NonAdmissible 


That {1, 2, 3,4,5,6,7,11,12,14,22,23,25} consists of forced gaps and that the set {0,18,19,20,21,24,26} 
can easily be confirmed by hand. 

Let us now check that 15 is non-admissible. If it was an element of a semigroup S € S{PF), then 
10(= 25-15) and 8(= 23-15) would be gaps of S. But then 17 is a big element, 13(= 23-10) is forced 
by exclusion (note that 25-10 = 15, 22-10 = 12 and 11-10 = 1 are gaps) and 9(= 23-14) is forced by 
exclusion too (25 - 14 = 11, 22 - 14 = 8 are gaps and 11 - 14 < 0). This is not possible, since 13 + 9 = 22 is 
a gap. Therefore 15 is non-admissible. 

5.3. A not so quick procedure to compute forced integers. Algorithm 2 is our procedure to compute 
forced integers that produces the best result in terms of the number of forced integers encountered. 
Besides Algorithm 1, it makes use of the function NonAdmissible. 

Remark 28. It is a consequence of Remark 21 that the second component of the list returned by Algo¬ 
rithm 2 is the set of small elements of a numerical semigroup. 


Input : PF 

Output: fail if some inconsistency is discovered; otherwise, returns [fg, fe ], where fg and fe are sets 
of forced gaps and forced elements, respectively, 
if type = 1 then 

^ return [Divisors(frob),[0,frob-i-l]]; 
fints := SimpleForcedIntegers(StartingForcedGaps(PF), [ ]); 
iifints = fail then 
return fail 
else 

if IsRange (Union(/z«ts)) then 
L return fints; 

nad := NonAdmissible(/z«is[l],/z«ts[2]); 
newgaps := Difference(«fld,/z«ts[l]); 

return Simp£eVocce6lntegecs(\]nion{newgaps, f ints[l]), f ints[2]); 

Algorithm 2: Forcedlntegers 


The correctness of Algorithm 2 follows from Proposition 26 and Theorem 20. 
Theorem 29. Algorithm 2 correctly produces the claimed output. 
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pseudo Frobenius numbers 

Time 

# F. Gaps 

# F. Elements 

QV 

NV 

QV 

NV 

QV 

NV 

[ 11, 22, 23, 25 ] 

2 

7 

13 

14 

7 

7 

[ 17, 27, 28, 29 ] 

2 

5 

16 

17 

8 

10 

[ 17, 19, 21, 25, 27 ] 

2 

9 

15 

16 

8 

8 

[ 15, 20, 27, 35 ] 

2 

10 

16 

16 

13 

13 

[ 12, 24, 25, 26, 28, 29 ] 

3 

11 

18 

22 

6 

9 

[ 145, 154, 205, 322, 376, 380 ] 

47 

1336 

82 

85 

52 

54 

[ 245, 281, 282, 292, 334, 373, 393, 424, 432, 454, 467 ] 

129 

2075 

116 

116 

53 

53 

[ 223, 434, 476, 513, 549, 728, 828, 838, 849, 953 ] 

213 

5866 

300 

318 

221 

253 

[ 219, 437, 600, 638, 683, 779, 801, 819, 880 ] 

205 

4838 

219 

224 

153 

161 

[ 103, no, 112, 137, 160, 178, 185 ] 

25 

262 

50 

51 

29 

31 


Table 1. Execution times for computing forced integers. 


5.4. Examples and execution times. Algorithm 1 can be used as a quick way to compute forced inte¬ 
gers. In fact, when called with the starting forced gaps and the empty set of elements, can be seen as a 
quick version of Algorithm 2. We use 

ForcedIntegers_OV(PF) 

as a shorthand for 

SimpleForcedIntegers(StartingForcedGaps (PF),[ ]) 

In the following example we use the names for the functions in the current implementation. 

Example 30. This example is meant to illustrate the difference between applying the quick and the nor¬ 
mal version of the algorithm. 

gap> pf := [ 103, 110, 112, 137, 160, 178, 185 ];; 
gap> fO := ForcedIntegersForPseudoFrobenius_OV(pf);;time; 

25 

gap> fN := ForcedIntegersForPseudoFrobenius(pf);;time; 

262 

gap> Length(f0[1]);Length(fN[1]);Length(fO[2]);Length(fN[2]); 

50 

51 
29 
31 

We have used the internal names in the package ForcedIntegersForPseudoFrobenius_OV and Forced- 
IntegersForPseudoFrobenius for ForcedIntegers_0V and Forcedlntegers, respectively. 

Table 1 collects some information concerning some execution times (as given by GAP) and the number 
of forced gaps and of forced elements both using Algorithm 1 (identified as QV (which stands for quick 
version)) and Algorithm 2 (identified as NV (which stands for normal version)). We observe that the 
execution times when using the quick version remain relatively small, even when the Frobenius number 
is large. 

Failure is usually detected very quickly, as should be clear and Table 2 somehow confirms. 

As one could expect, there are examples where failure is not detected with the quick version. 
gap> pf := [ 25, 29, 33, 35, 38, 41, 46 ]; 

[ 25, 29, 33, 35, 38, 41, 46 ] 

gap> ForcedIntegersForPseudoFrobenius(pf); 

fail 

gap> ForcedIntegersForPseudoFrobenius_OV(pf); 

[ [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 16, 17, 19, 21, 23, 25, 29, 33, 35, 
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pseudo Frobenius numbers 

Time 

QV 

NV 

[ 18, 42, 58, 88 , 94 ] 

6 

6 

[ 20, 27, 34, 35, 37, 42, 48, 80 ] 

2 

3 

[ 30, 104, 118, 147, 197, 292, 298, 315, 333, 384, 408 ] 

32 

43 

[ 36, 37, 219, 233, 304, 410, 413, 431, 438, 458 ] 

35 

32 

[ 89, 411, 446, 502, 557, 600, 605, 631, 636, 796, 801, 915 ] 

223 

233 

[ 56, 134, 136, 137, 158, 248, 277, 373, 383, 389, 487, 558, 566, 621, 691, 825, 836 ] 

103 

113 


Table 2. When failure occurs... 


38, 41, 46 ], [ 0, 30, 34, 36, 37, 39, 40, 42, 43, 44, 45, 47 ] ] 
gap> pf := [ 22, 23, 24, 25, 26 ]; 

[ 22, 23, 24, 25, 26 ] 

gap> ForcedIntegersForPseudoFrobenius_OV(pf); 

[ [ 1, 2, 3, 4, 5, 6, 8, 11, 12, 13, 22, 23, 24, 25, 26 ], [ 0, 20, 21, 27 ] ] 

gap> ForcedIntegersForPseudoFrobenius(pf); 

fail 

On the other hand, we have not been able to detect any set PF candidate to be the set of pseudo- 
Frobenius numbers of a numerical semigroup that passes the normal version and such that 5(PF) = 0. 
Despite the various millions of tests made, we do not feel comfortable on leaving it as a conjecture; we 
leave it as a question instead. 

Question 31. If Algorithm 2 with input PF does not return/ad, then 5(PF) ^ 0? 

We observe that in view of the execution times illustrated in Table 1, a positive answer to Question 31 
would imply that Algorithm 2 constitutes a satisfactory answer to Question 1. 

6. An approach based on forced integers 

In this section we present our main algorithm (Algorithm 3), which computes 5(PF). Its correctness is 
stated in Theorem 39, whose proof is built from almost all the results preceding it in the paper. 

After considering some initial cases, the algorithm makes a call to RecursiveDepthFirstExploreTree 
which is a recursive function used to construct a tree whose nodes are labeled by pairs (X,Y) where X 
is a list of forced gaps and Y is a list of forced elements. Thus we implicitly have lists of free integers in 
each node: the complement of X U Y in the set U = {!,...,g„}, where PF = {gi < ••• < g„}. Nodes with an 
empty set of free integers are the leafs in our tree. 

A node (X, Y) such that there exists a numerical semigroup S € 5(PF) for which X c gaps(S) and 
Y c smalls(S) is said to be PV-admissible, or simply admissible, if PF is understood. A node that is not 
admissible is called a dead node. 

Remark 32. The knowledge of some forced integers allows us to identify immediately some dead nodes: 
if (G,E) is such that G consists of forced gaps and E consists of forced elements, then any node (X, Y) 
such that Xn£ 7 t 0 orYnG 7 t 0 isa dead node. 

Remark 33. Let (X, Y) be a leaf that is not a dead node. It follows from the construction (see Remarks 21 
and 28) that there exists a numerical semigroup S such that (gaps(S),smalls(S)) = (X, Y). 

Remark 34. Observe that if PF = or PF = {gi/2 < gi), then the set of numerical semigroups with 
pseudo-Frobenius numbers PF corresponds with the set of irreducible numerical semigroups having 
Frobenius number gi, see Appendix A. In this case we will use the fast procedure presented in [1]. 

6.1. The recursive construction of a tree. A naive idea is to start with a list of free integers (for some 
set of forced integers) and turn each one of these free integers into either a gap or an element. Assuming 
that the number of free integers is n, the number of possibilities is 2 ", thus checking each of these 
possibilities for being in correspondence with the set of gaps of a numerical semigroup with PF as set of 
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pseudo-Frobenius numbers is unfeasible, unless n is small. Nevertheless, this naive idea can be adapted 
and take advantage of the already observed facts that elements can force gaps and vice-versa. Although 
there are examples for which fixing one integer does not force anything else, which let us expect that 
nothing good can result from a worst case complexity analysis, in practice it works quite well. We give 
some examples, but leave a detailed analysis of the complexity (perhaps average complexity) as an open 
problem. 

The procedure we use is, for each integer v in the current list of free integers, compute all numerical 
semigroups containing v and the forced elements, and having the forced gaps as gaps. We proceed 
recursively and use backtracking when a semigroup or a dead node is found. When we revisit the node, 
we then suppose that v is a gap and continue with the next free integer. 

Before proceeding with the pseudo-code for the recursive function RecursiveDepthFirstExploreTree 
that constructs the tree in a depth first manner, let us give an example which we accompany by a picture 
(Figure 2). 



Figure 2. The numerical semigroups with pseudo-Frobenius numbers {15,20,27,35}. 

Example 35. Let PF = {15,20,27,35}. From Example 24 we have a pair of lists of forced gaps and forced 
integers, which leaves the list 

(2) f = {11,13,14,17,18,21,22,24} 

of free integers. 

The leaves contained in the branch that descends from 11 in Eigure 2 consist of the semigroups con¬ 
taining 11 as forced element. 

All the remaining semigroups in 5(PF) must have 11 as gap. The process then continues as Figure 2 
illustrates: considering 13 as a forced integer, developing its corresponding subtree and so on. 

Let us now look at the integer 21 in the root of the tree. At this point, all the semigroups in 5(PF) 
containing some integer preceding 21 in f have been computed. Thus, any new semigroup must have 
the elements in {11,13,14,17,18} as gaps. From the computations one concludes that no new semigroup 
appears having 21 as an element. Thus, if some not already computed semigroup exists fulfilling the 
current conditions, then it must have 21 as gap. One can check that this can not happen. The remaining 
elements (22 and 24) need not to be checked, since 21 had to be either an element or a gap. Therefore 
we can stop. 
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Next we give pseudo-code for RecursiveDepthFirstExploreTree. The recursion ends when no more 
than one free integer is left. A call to the function EndingCondition is then made. 

The variable semigroups works as a container which serves to store all the semigroups we are looking 
for, as they are discovered. 


EndingCondition(g,e) # g and e are such that #Lf\(gUe)<l 

# g and e represent lists of gaps and elements, respectively 
free := ShallowCopy(Difference([l..frob],Union(g,e))); 

1 if Length(/rec) = 0 then 

2 if First(Dif f erence(g, PF),p/ ^ Intersection(p/ -i- Dif f erence(e, [0]), g) = 0) = fail then 
^ Add(semigroups, NumericalSemigroupByGaps(g)); 

return; 


3 

4 


if Fength(/ree) = 1 then 

if RepresentsGapsOfNumericalSemigroup(g) then 

if First(Difference(g, PF),p/ ^ Intersection(p/-i-Dif f erence(Union(e,/ree), [0]),g) = 
0) = fail then 

1^ Add(semigroups, NumericalSemigroupByGaps(g)); 


if RepresentsGapsOfNumericalSemigroup(Union(g,/ree)) then 
5 if First(Diff erence(g, PF),p/ ^ Intersection(p/ -i- Dif f erence(e, [0]), Union{g, f ree)) = 

0) = fail then 

|_ Add(semigroups,NumericalSemigroupByGaps(Union(g,/ree))); 

return; 


Function EndingCondition 


Lemma 36. Function EndingCondition either does nothing or adds to semigroups a numerical semigroup S 
such that PE(S) = PF. 

Proof. It suffices to observe that any of the conditions in the "if" statements of Lines 2, 4, 5 guarantee 
that no forced gap outside PF can be a pseudo-Erobenius number. □ 

Notice that when EndingCondition does nothing, it is because one of the following reasons: 

• there is no numerical semigroup whose set of gaps is the first component of the input, 

• the resulting semigroup does not have PF as set of pseudo-Erobenius numbers (it actually has 
more pseudo-Erobenius numbers), 

• there is a free element that cannot be neither a gap nor an element. 

Observe that recursion in RecursiveDepthEirstExploreTree ends (in Line 3) when there is at most one 
free element. 

Proposition 37. Let {fg,fe) be a pair of disjoint sets of integers contained in U = {0,...,g„ + 1}. After 
the execution of the function RecursiveDepthFirstExploreTree with input (fg,fe), semigroups contains all 
numerical semigroups S such that S € 5(PF), fg c gaps(S) and fe c smalls(S). 

Proof Denote by A the set of all numerical semigroups S such that fg c gaps(S), fe c smalls(S). We 
have to prove that semigroups n A = 5(PF) n A. 

c. It suffices to observe that the numerical semigroups are added to semigroups by the function End¬ 
ingCondition and these belong to 5(PF), by Lemma 36. 

3. Let S € 5(PF) be such that fg c gaps(S) and fe c smalls(S). If #U\ (fgCfe) e {0,1}, then the function 
EndingCondition is called and it gives the correct output. Otherwise, we enter a recursion. 

We will prove by induction on the number of free elements that the output is correct. Let us consider 
the smallest integer v eU \ {fg U fe). RecursiveDepthEirstExploreTree is called with input {fg,fe), and 
will enter the while loop, adding v to fe and computing new forced integers (SimpleForcedlntegers). 
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RecursiveDepthFirstExploreTree([fg,fe]) # fg and fe consist of forced gaps and 

# forced elements, respectively 
current free := ShallowCopy(Dif f erence([l..frob], Union(fg,fe))); 

«/g := ShallowCopy(fg); # used to store new forced gaps 

# without creating conflicts in the memory 

while Length(curre«f/rer) > 1 do 
V := currentfree[\\, 

1 left := SimpleForcedIntegers(«/g,Union(fe, [v])); 

2 if/e/t = fail then 

right := SimpleForcedIntegers(Union(«/g, [v]),fe); 
ii {right = fa\[) or {lntecsection{right[l],right[2]) then 
L break; 

else 

^ RecursiveDepthFirstExploreTree {[left[l]Jeft[2]]); 
nfg := Union(«/g, [v]); 

current free := Diiierence{currentfree, [v]); 

3 if Length{currentfree) < 1 then 
^ EndingCondition (fg,fe); 

Function RecursiveDepthFirstExploreTree 

• If V € S, then left in Line 1 will not be equal to fail (Theorem 20), and we will call Recur¬ 
siveDepthFirstExploreTree with a lager set fe, having in consequence less free integers. By in¬ 
duction, S is added to semigroups. 

• Now assume that v € S. After the execution of the if-then-else starting in Line 2, v is added to 
the set of gaps. We have then one element less in the list of free integers, fg U {v} c gaps(S) 
and fe c smalls(S), whence the effect is the same as if we called RecursiveDepthFirstExploreTree 
with arguments fg U {v} c gaps(S) and fe c smalls(S). By induction hypothesis, S is added to 

semigroups. □ 

Observe that if 5(PE) ^ 0 and the sets fg and fe considered in Proposition 37 consist of forced gaps 
and forced elements, respectively, then A c 5(PE). Therefore we have proved the following corollary. 

Corollary 38. If the function RecursiveDepthFirstExploreTree is called with a pair {fg,fe), where fg consists 
of forced gaps and fe consists of forced elements, then at the end of its execution we have that semigroups = 
5(PE). 


Input : PE 
Output: the set 5(PE). 

1 if (type = 1) or (type = 2 and PE[1] = PE[2]/2) then 

^ return IrreducibleNumericalSemigroupsWithFrobeniusNumber(frob); 

2 if not (PE[type] - PE[type - 1] > PE[1]) then 
L return 0; 

3 root:=EorcedIntegers (PE); 
if root = fail then 

L return 0; 
semigroups := 0; 

RecursiveDepthE irstExploreTree(root); 

return semigroups; 


Algorithm 3: NumericalSemigroupsWithPseudoFrobeniusNumbers 
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6.2. The main algorithm. We observe that Algorithm 3 is not efficient when there are many free integers 
(for some set of forced integers). A possible attempt to improve it could be to replace the recursive 
function by a much more efficient depth first search algorithm to parse the tree in question. Another 
one is to use other sophisticated theoretical tools. What is done in Appendix A could be seen as an 
attempt. Having at our disposal more than one approach, we can take advantage of choosing the most 
efficient for each situation. 

Theorem 39. Let PF = < ••• < g„} be a set of positive integers. The output of Algorithm 3 with input PF is 

5(PF). 

Proof Line 1 is justified by Remark 34. The justification for Line 2 comes from Corollary 10. If the stated 
necessary condition is not fulfilled, then 5(PF) = 0, which is precisely the set returned by the algorithm. 
As Algorithm 2 returns/flz7 precisely when it is found some forced gap that is at the same time a forced 
element. Proposition 3 assures us that 5(PF) = 0, which is the set returned by the algorithm. 

When nothing of the above holds, the variable semigroups is initialized as the empty set and it is made 
a call to the recursive function RecursiveDepthPirstExploreTree. As we are considering this variable 
global to the functions EndingCondition, the result follows from Corollary 38. □ 

7. Running times and examples 

The number of semigroups can be quite large compared to the number of free elements. The following 
example illustrate this possibility. 

Example AO. gap> pf := [ 68, 71, 163, 196 ];; 
gap> forced := ForcedIntegersForPseudoFrobenius(pf);; 
gap> free ;= Difference([1..Maximum(pf)],Union(forced));; 
gap> Length(free); 

38 

gap> list ;= NumericalSemigroupsWithPseudoFrobeniusNumbers(pf);; 
gap> Length(list); 

1608 

In the continuation of previous GAP session we do a kind of verification of the result. 
gap> ns := RandomList(list); 

<Numerical semigroup> 

gap> MinimalGeneratingSystem(ns); 

[ 35, 38, 65, 81, 89, 94, 99, 101, 104, 106, 109, 110, 112, 113, 117, 118, 

121, 122, 133 ] 

gap> PseudoFrobeniusOfNumericalSemigroup(ns); 

[ 68, 71, 163, 196 ] 

Table 3 is meant to illustrate some timings. Its content was produced by using repeatedly the com¬ 
mands of the first part of the previous example. The candidates to PF were obtained randomly. We 
observe that, although depending on some factors (such as the type or the Erobenius number we are 
interested in), the vast majority of the candidates would lead to the empty set. We do not consider them 
in this table (some were given in Table 2). 


8. Random 

Sometimes one may just be interested in obtaining one numerical semigroup with PF as set of pseudo- 
Erobenius numbers. Algorithm 3 may be too slow (it gives much information that will not be used). One 
could adapt the algorithm to stop once it encounters the first semigroup, but the information had to be 
transmitted recursively and one would end up with a slow algorithm. Next we propose an alternative 
(Algorithm 4), the first part of which is similar to the initial part of Algorithm 3. The main difference is 
in the usage of the function AnIrreducibleNumericalSemigroupsWithErobeniusNumber instead of the 
function IrreducibleNumericalSemigroupsWithErobeniusNumber, both available in the numericalsgps 
package. Eor the second part, instead of calling the recursive function, it tries to gess a path that leads 
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pseudo Frobenius numbers 

# free 

# semigroups 

time 

[ 15, 27, 31, 43, 47 ] 

0 

1 

3 

[ 16, 30, 33, 37 ] 

9 

3 

12 

[ 40, 65, 80, 89, 107, 110, 130 ] 

5 

3 

29 

[ 32, 35, 44, 45, 48 ] 

13 

7 

24 

[ 40, 65, 89, 91, 100, 106 ] 

24 

9 

99 

[ 36, 50, 56, 57, 63 ] 

25 

39 

123 

[ 43, 50, 52, 65 ] 

35 

213 

605 

[ 68, 71, 163, 196 ] 

38 

1608 

16603 

[ 38, 57, 67, 74, 79 ] 

40 

155 

527 

[ 68, 72, 76, 77 ] 

46 

177 

607 

[ 62, 78, 99, 129, 130 ] 

53 

4077 

28622 

[ 128, 131, 146, 151, 180, 216, 224, 267, 271, 287 ] 

54 

954 

24253 

[ 84, 103, 144, 202, 230, 242, 245 ] 

56 

14292 

277094 

[ 66, 85, 86, 92 ] 

55 

950 

4683 

[ 76, 79, 88, 102 ] 

64 

1409 

6505 

[ 61, 67, 94, 105 ] 

69 

4432 

21471 

[ 114, 150, 179, 182, 231, 236, 254, 321 ] 

69 

302929 

7121020 

[ 62, 73, 166, 190, 203 ] 

77 

9934 

134554 

[ 102, 104, 118, 123, 134, 146, 149 ] 

87 

15910 

149910 


Table 3. Some examples of execution data of the main algorithm. 


to a leaf. Starts choosing at random a free integer v and tests its non admissibility (by checking whether 
SimpleForcedlntegers returns fail when called with v as if it was forced). If one does not conclude 
that V is non admissible, it is assumed to be a forced integer. There is an option that is part of the 
implementation to give a bound for the maximum number of attempts the function does. Its usage is 
ilustrated in Examples 41 and 42. 

It may happen that no semigroup has the given set as set of pseudo-Frobenius elements, and thus the 
output will simply be fail. 

Example 41. We look for a random numerical semigroup with PF = {100,453,537,543}. The first execu¬ 
tion of the function yields: 
gap> pf := [ 100, 453, 537, 543 ];; 

gap> ns := RandomNumericalSemigroupWithPseudoFrobeniusNumbers(pf);;time; 

MinimalGeneratingSystem(ns); 

2440 

[ 66 , 94, 106, 126, 166, 184, 194, 206, 209, 216, 224, 230, 235, 246, 256, 263, 

267, 284, 295, 309, 363, 374, 379, 385, 391, 413 ] 

While if we execute it a second time we get another semigroup with the desired set of pseudo-Frobenius 
numbers (clearly, the reader might obtain different outputs): 
gap> pf := [ 100, 453, 537, 543 ];; 

gap> ns := RandomNumericalSemigroupWithPseudoFrobeniusNumbers(pf);;time; 
MinimalGeneratingSystem(ns); 

7302 

[ 94, 106, 123, 134, 162, 178, 184, 194, 204, 206, 222, 223, 234, 235, 238, 248, 

251, 262, 263, 266, 270, 276, 279, 283, 293, 304, 313, 336, 348, 367, 383, 415 ] 

Example 42. If one of the free integers can neither be a gap nor an element, no semigroup exists. 
gap> pf := [ 30, 104, 118, 147, 197, 292, 298, 315, 333, 384, 408 ];; 
gap> ns := RandomNumericalSemigroupWithPseudoFrobeniusNumbers( 

> rec(pseudo_frobenius := pf, maxattempts := 100));time; 
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Input : PF, max.attempts 

Output: One numerical semigroup S (at random) such that PF(S) = PF if it discovers some, fail if it 
discovers that no semigroup exists... Otherwise, suggests to use more attempts 
if (type = 1) or (type = 2 and PF[1] = PF[2]/2) then 
^ return AnlrreducibleNumericalSemigroupsWithFrobeniusNumber(frob); 

if not PF[type] - PF[type - 1] > PF[1]) then 
L return fail; 

fjnts := Forcedlntegers(PF); 
iifjnts = fail then 
L return fail; 

free := Dif f erence([l..frob], Union(/Juts)); 
for i € [l..max-attempts] do 
while free 0 do 

V := RandomList(/ree); 

nfig := SimpleForcedIntegers(Union(/J«ts[l], [v]), f Jnts[2]); 
nfie := SimpleForcedIntegers(/J«ts[l], Union)/J«ts[2], [v])); 
if nfig ^ fail then 

if IsRange(Union(«/zg)) then 
I return NumericalSemigroupByGaps(«/zg[l]); 
else 

|_ free := Difference([l..frob],Union(«/z^)); 

else 

if nf ie ^ fail then 

if IsRange(Union(«/ze)) then 
I return NumericalSemigroupByGaps(«/ze[l]); 

else 

/ free := Difference([l..frob],Union(«/ze)); 

else 

1_ break; 

# Info: Increase the number of attempts... 

Algorithm 4: RandomNumericalSemigroupWithPseudoFrobeniusNumbers 


fail 
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Appendix A. An approach based on irreducible numerical semigroups 

We present here an alternative way to compute the set of numerical semigroups with a given set of 
pseudo-Frobenius numbers. In general, this procedure is slower than the presented above, though we 
have not been able to characterize when this happens. We include it the manuscript since it was the 
initial implementation and was used to test the other one. 

A numerical semigroup is irreducible if it cannot be expressed as the intersection of two numerical 
semigroups properly containing it. It turns out that a numerical semigroup S is irreducible if and only 
if either PF(S) = {F(S)} or PF(S) = {F(S)/2,F(S)} (see [ 8 , Chapter 3]). Irreducible numerical semigroups 
can be also characterized as those maximal (with respect to set inclusion) numerical semigroups in the 
set of all numerical semigroups with given Frobenius number. 

The maximality of irreducible numerical semigroups in the set of all numerical semigroups with given 
Frobenius number implies that every numerical semigroup is contained in an irreducible numerical with 
its same Frobenius number. Actually, we can say more. 

Lemma 43. Let S be a numerical semigroup. There exists an irreducible numerical semigroup T such that 




(1) F(S) = F(T), 

(2) ]F(S)/2,F(S)[nPF(S)cT. 
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Proof. Let / = F(S) and let F =]//2,/[nPF(S). We claim that S' = S U f is a numerical semigroup with 
F(S') = /. Take s,s'€S'\{0}. 

• If both s and s' are in S, then s + s' € S c S'. 

• If s € S and s' € F, then s + s' € S, because s' € PF(S). 

• If s,s' € F, then s + s' > /, and so s + s' € S c S'. 

Let us show that F(S') = /. Assume to the contrary that this is not the case, and consequently / € S'. 
Then as all the elements in F are greater than f/2, and f €S, there must be an element s € S and g € PF(S) 
such that g + s = f. But this is impossible, since all elements in PF(S) are incomparable with respect to 
<S (Lemma 4). 

If S' is not irreducible, then as irreducible numerical semigroups are maximal in the set of numerical 
semigroups with fixed Frobenius number, there exists T with F(S) = F(S') = F(T) and containing S'; 
whence fulfilling the desired conditions. □ 

With all these ingredients we get the following result. 

Proposition 44. Let S be a nonirreducible numerical semigroup with PF(S) = {gi < ■■■ < gj^^], k > 2. Then 
there exists a chain 

S = So C Si = S U {xi} C ••• C S/ = S U {xi,...,x/}, 

with 

(1) Si irreducible, 

(2) Xi = max(S; \ S,_i)/or all i e {1,...,/}, 

(31 (gk/2,gk)nPF(S)cSi, 

(4) for every z € {1,...,/}, x,- is a minimal generator of Si such that gj - x, € Sjfor some j 

(5) for every z € {1,..., /} and f e PF(S;) with f + gi., there exists j e {I,...,k-l} such that gj- f ^ Si- 

Proof Let T be as in Lemma 43. Construct a chain joining S and T by setting Sq = S and S; = S,_i U {x,}, 
with X; = max(T \ S;_i). Then S, is a numerical semigroup, and X; is a minimal generator of S, and a 
pseudo-Frobenius number for S,_i ([ 8 , Lemma 4.35]). Since the complement T \ S is finite, for some k, 
Sk = T. 

Clearly, Xi = gj., and thus Xi -gjt = 0 ^ S. Now let z € {2,...,k}. Then x,- € T \ S, and by Lemma 4, there 
exists j e {I,...,k} such that gj -x,- € S c S,-. 

Take / € PF(S,) \ {g^}- Then / g S and thus gj - f e S for some j e {l,...,k} (Lemma 4 once more). 
Consequently gj- f ^ S,. Notice that j + k, since f -gk<^- D 

Given a candidate set PF of pseudo-Frobenius numbers with maximum element /, we can use the 
above procedure to construct from the set of all irreducible numerical semigroups with Frobenius num¬ 
ber /, the set of all numerical semigroups having PF as a set its pseudo-Frobenius numbers. In order to 
compute the set of all irreducible numerical semigroups with Frobenius number / we use implementa¬ 
tion of the procedure presented in [1] that is already part of [4]. We have slightly modified the algorithm 
in [ 1 ] to compute the set of irreducible numerical semigroups containing a given set of integers, and 
these integers are the first component of Forcedlntegers (PF). For every irreducible element in the list 
we then remove those minimal generators fulfilling condition (4) in Proposition 44. We add to our list 
of semigroups the semigroups obtained in the preceding step for which condition (5) holds, and then we 
proceed recursively. 

Example 45. Let us illustrate the above procedure with PF = {10,13}. The number of irreducible numer¬ 
ical semigroup with Frobenius number 13 is 8 . However, if we first call Forcedlntegers we get: 

gap> ForcedIntegersForPseudoFrobenius([10,13]); 

[ [ 1, 2, 3, 5, 6, 10, 13 ], [ 0, 7, 8 , 11, 12, 14 ] ] 
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Since we have modified the function IrreducibleNumericalSemigroupsWithFrobeniusNumber to 
output only those irreducible numerical semigroups containing {0,7,8,11,12,14}, we obtain only two 
irreducible numerical semigroups: Si = (4,7,10) and S 2 = (7,8,9,10,11,12). 

For Si the only minimal generator that fulfills the conditions in Proposition 44 is 10. If we remove 10 
from Si, we obtain Ti = (4,7,17), which already has the desired set of pseudo-Frobenius number. 

As for S 2 , again 10 is the only minimal generator fulfilling the conditions in Proposition 44, and we 
obtain T 2 = (7,8,9,11,12). This semigroup has pseudo-Frobenius number set equal to PF, and so, as with 
Ti we do not need to look for new minimal generators to remove. 

Thus, the only numerical semigroups with pseudo-Frobenius number set {10,13} are Ti and T 2 . 

References 

(1] V. Blanco, J. C. Rosales, The tree of irreducible numerical semigroups with fixed Frobenius number. Forum Math. 25 
(2013), 1249-1261. 

(2] H. Bresinsky, Symmetric semigroups of integers generated by 4 elements, Manuscripta Math. 17 (1975), 205-219. 

(3] M. Delgado, "intpic", a CAP package for drawing integers. Available via http: / /www .gap-system, org/. 

{4} M. Delgado, R A. Garcia-Sanchez and J. Morals, "NumericalSgps", a GAP package for numerical semigroups, Version 0.990; 
2015. Available via http: / / www. gap -system.org/. 

(5] The GAP Group, GAP - Groups, Algorithms, and Programming, Version 4.7.7; 2015. Available via http://www. 
gap-system. org/. 

(6] J. Komeda, On the existence of Weierstrass points with a certain semigroup generated by 4 elements, Tsukuba J. Math. 6 
(1982), 237-270. 

{7} A. M. Robles-Perez and J. C. Rosales, The genus, the Frobenius number, and the pseudo-Frobenius numbers of numerical 
semigroups with type two, preprint. 

(8] J. C. Rosales and P. A. Garcia-Sanchez, "Numerical Semigroups", Developments in Maths. 20, Springer (2010). 

(9] J. C. Rosales, P. A. Garcia-Sanchez, J. 1. Garcia-Garcia, and J. A. Jimenez-Madrid, Fundamental gaps in numerical semi¬ 
groups, J. Pure Appl. Algebra 189 (2004), 301-313. 

CMUP, Departamento de Matematica, Faculdade de Ciencias, Universidade do Porto, Rua do Campo Alegre 687, 4169- 
007 Porto, Portugal 

E-mail address: mdelgacloQf c .up.pt 

Departamento de Algebra, Universidad de Granada, 18071 Granada, Espana 
E-mail address: pedro@ugr.es 

Departamento de Matematica Aplicada, Universidad de Granada, 18071 Granada, Espana 
E-mail address: aroblesQugr. es 



